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We study the combined counting statistics of two capacitively coupled transport channels. In 
particular, we examine the conditions necessary for utilizing one channel as detector sensitive to the 
occupation of the other. A good detector fidelity may be achieved in a bistable regime when the 
tunneling rates through the two channels are vastly different - even when the physical back-action 
of the detector on the probed channel is large. Our methods allow to estimate the error of charge 
counting detectors from time-resolved current measurements - which have been obtained in recent 
experiments - alone. 



Quantum transport is a generic example of nonequi- 
librium quantum dynamics^. A low-dimensional quan- 
tum system (e.g. quantum dots^, molecules^, nanotubes^ 
etc.) is typically coupled via particle exchange to mul- 
tiple reservoirs held at different equilibria. The detec- 
tion of stochastic particle transfer into a reservoir yields 
the full counting statistics (FCS)£. The FCS provides 
a tool to access system properties via an indirect mea- 
surement. The technique of n- resolved master equations^ 
conveniently allows to extract the FCS from microscopic 
models (see e.g<£). Recently, the value of this approach 
has been demonstrated experimentally by capacitively 
coupling a single quantum dot (QD) to a quantum point 
contact (QPC)£. Here, we demonstrate that in principle, 
also two capacitively coupled QDs should show similar 
behavior as the QPC-QD configuration in certain param- 
eter limits. However, the simplicity of our model enables 
us to estimate the physical detector back-action and the 
detection error from the FCS. 

We consider a system composed of two nearby two- 
terminal single resonant level systems^ A and B (com- 
pare Fig. [Th,) H s = e A d\d A + e B d B d B + Ud\dAd B d B , 
where &a/b denote the level energies of the sin- 
gle levels and U models Coulomb interaction. The 
system is coupled to four fcrmionic reservoirs via 



H SB = Efca tka,A d Ac\ aA + J2ka tka,Bd B c\^ B + h.C, 

where t^A/B denote tunneling rates to the adjacent 
leads a £ {L, R}. We will use subsystem A as a detector 
for the state of subsystem B. Transport through each 
level is enabled by applying nonvanishing bias voltages, 
and the two transport channels influence each other by 
the Coulomb interaction U. 

Since we are interested in the charge FCS not only at 
one junction but at the interplay of the two channels, we 
introduce two virtual detectors in the tunneling terms as- 
sociated with right leads (see also Ref>^). Following the 
Born- Markov-Secular— approximation scheme, we arrive 
at an (n A , JT-s)-resolved master equation of the form 
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where = (n A , n B \ p \ua, n B ), which couples 

only the diagonals of the density matrix to each 



other. The total system (QDs and virtual detec- 
tors) density matrix can at all times be written as 

P(t) = J2n A , nB P (nAnB) (t) ® \n A )(n A \ ® \n B )(n B \, 
such that the probability to measure ha tunneled 
particles in the detector channel and n B tunneled 
particles in the system channel after time t is given by 
Pn A n B {t) = Tr {p {nAnB \t)} (measurement postulate). 
Performing a two-dimensional Fourier-transform of 



Eq. © Via p(x,£,t) = En A ,n B P [nAnB) 

we obtain a Fourier-transformed 
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Liouvillian 

P(X)£:*) = £(X)£)p(X)£j*) W1 th two counting fields 
C( X ,0 = £oo + e +iX £+o + e- l *£-o + e+ l «A,+ + e^A>-. 
This equation can be solved formally, and for a given 

initial condition p(0) = p 00) <8 |0) (0| ® |0) (0|, the density 
matrix p( nAnB \At) can be obtained from the inverse 
Fourier transform. With the "superjump" superoperator 
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FIG. 1: (Color Online) (a) Spatial sketch of the physical 
setup. Tunneling is described by rates F and 77, respectively. 
Detectors are placed at the right junctions - with counting 
fields x an d £- Coulomb interaction U leads to an effective 
shift of the dot levels (upper orange dots in panels b-f ) , which 
influences the current through each channel, (b-f) Band- 
scheme sketches with the associated Fermi functions f^j^(to) 
(labels in panel b) for different bias and temperature config- 
urations (see text). 



2 



we obtain after time t for the density matrix 
p (nA,n B )(£) = j(nA,n B )^ p (oo)^ Measuring tnc number 

of tunneled particles through both channels at this point 
will collapse the density matrix again. The explicit form 
of the Liouvillian considered in this work is given by a 
special bistable case of a size-scalable variant - with an 
additional counting field £. 

We consider only the case of unidirectional trans- 
port as shown in Figs. [T|3-f: All right-associated Fermi 
functions vanish at the energy scales of the system 
f R ( €j ) = f R (ej +U) = (with j G {A, £?}), and the Li- 
ouvillian is given by = TCa(x) + V^-b(£,) with 
the flat tunneling rates r = 2ir \tk a ,A\ 2 3(u — ^>k.a) 
and j] = 2n^ k \tka,s\ 2 S(uj — uik,a)- In the basis 
jOoo,oo,Pio,io,Poi,oi,Pii,n, the superoperators read 
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■U. Bistability 

occurs when e.g. t] —> and (ca) ^ /l (ca), since the 
block structure of the Liouvillian supports two different 
currents (and stationary states) in this case^. 

Current Trajectories. When channel B is at infi- 
nite bias and at sufficiently large temperatures k^T > U, 
the left-associated Fermi functions f£ (ca) and /^(ea) 
will assume intermediate values between zero and one 
(see Fig [I}}), which yields two different currents through 
A depending on whether B is occupied or not. Ignor- 
ing the number of transferred charges through channel 
B, the probability of obtaining ha charges after time 
At equates to P„ A (At) = Tr{£„ B J^ n ^(At)p } = 

Trj^/^^^A^po}- That is, to obtain a trajectory of cur- 
rent measurements performed equidistantly at intervals 
At numerically, a random number according to the dis- 
tribution P nA (At) must be drawn. The outcome of this 
then corresponds to a measurement of ua particles (with 
current I a = n A /At). For simplicity, we exploit the 
translational invariance of Eq. ([1]) and shift ua to zero 
after each measurement. However, for the subsequent 
evolution, one now has to use the normalized density 
matrix as initial condition, which leads to a tem- 

poral correlation of the measured currents, see Fig. [2] In 
principle, also for unidirectional transport the number of 
different outcomes SEJl G {0, 1, 2, . . .} is infinite, but fortu- 
nately, there exists a natural cutoff as for ua > TAt the 
operators J^ lA \At) become exponentially small. The 
trajectory in Fig. [5] is very similar to results obtained 
for QPCs^. The blips from low to high currents (ar- 
rows) should therefore indicate single charges tunneling 
from channel B to its right junction. However, as tun- 
neling through channel A is also a stochastic process, 
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FIG. 2: (Color Online) Trajectory for configuration in Fig. [Tb 
for 1000 current measurements at time intervals of At. Back- 
ground lines display the currents and widths (shaded regions) 
for the single resonant level evaluated at the effective QD level 
one obtains for empty (high current) or filled (low current) 
channel B. Temperature and chemical potentials chosen such 
that /l(£a) = 0.612 and /l(<Ta) = 0.181. Inset (parameters 
TAt = 10, rjAt = 0.1, At = 100.0 arb. u.) shows the corre- 
lation of both (coarse-grained) currents for the configuration 
in Fig. When the detector current (solid black) rises, one 
most often also observes a spike (single charge) in the cur- 
rent through channel B (dashed red, rescaled). However, due 
to stochastic behavior of the detector, jump events may be 
missed completely (bold red arrow). 



jump events may in principle be missed (compare also 
the inset). Therefore, this raises the question of detector 
fidelity and back-action. In the following, we will inves- 
tigate this for configurations (c)-(f) in Fig.Q] 

Unperturbed Level (Fig. \J]c). In this limit, the 
current through A is blocked completely, such that 
the detector is turned off. With the replacements 
= fi(e A ) = and f£(e B ) = /£(e B ) = 1 in 
Eq. ([3]), we obtain for the Laplace transform M(x, £, z) = 
Tr — C(XiO] 1 p| °f the moment-generating func- 
tion (MGF) 



2[z 2 +2rjz- (e« - l)rf 



(4) 



from which one may analytically^ deduce the probability 
distribution of a single resonant level in the symmetric 
infinite bias limit. 

Infinite Bias (Fig. \Jjfl). When both channels are 
held at infinite bias - replacements (ca) = ft (^a) = 1 
and /f (e B ) = /f (e B ) = 1 in Eq. © - the detector 
cannot distinguish the different states of the probed sys- 
tem. From the exact matrix exponential of the Liou- 
villian we obtain P nA n B (t) = Pn A (t)P nB (t), where each 
probability distribution corresponds to that of a single 
resonant levels (with r\ — > T for channel A). In Laplace 
space, this becomes visible by considering the reduc- 
tion to Eq. @ via M(0, £, z) = M B (£, z) and similarly 
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M(x, 0, z) = M A (x,z). Clearly, we obtain for the first 
and second (long-term limit) cumulants of the detector 
(current and noise) {{n A }) = an d((n^)) —> ■ 
for the probed system ((ng)) = \ and ((n 2 B )) 
From the factorization of the probability distributions it 
follows also that the cross-correlations vanish at all times 
((n A n B )) = (n A n B ) ~ (n A ) (n B ) = 0. 

Charge Detector Limit (Fig.[I]p). When f£(eA) = 
1, ft (ca) = 0, and ff(e B ) = /£ (eg) = 1 in Eq. © - the 
detector is bistable^, which can be used to distinguish 
the two different states of channel B: When channel B 
is occupied, current through A is blocked completely, 
whereas it is maximal (at its infinite bias value) other- 
wise. We find that the Laplace transform of the MGF 
fulfills Af(0,£,z) = M B (£,z) - compare Eq. Q - which 
demonstrates that here the physical action of the detector 
on the probed system vanishes. In contrast, the physi- 
cal influence of the system on the detector is not negli- 
gible, which is demonstrated by M(x,0,z) 7^ Ma(Xi z )- 
Therefore, the role of the detector will be purely informa- 
tional (measurement postulate) without physical back- 
action on the system. 

We obtain for the first cumulants ((n A )) = 
((rtg)) = whereas the second cumulants asymptot- 
ically approach the long-term limits 
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i.e., there exists a slight negative cross-correlation be- 
tween the two channels, which surprisingly saturates. 

In order to operate the channel A as a detector for 
the FCS of channel B, we need to scan channel A at 
much larger rates than B (T 3> 77). In particular, to 
obtain a meaningful detector current [associated with a 
transiently bimodali 2 - distribution P nA (At)], we require 
TAr; 3> 1 and to see the slow switching in time-resolved 
current measurements it is also necessary that 77 At <C 1. 
Then, measuring a vanishing detector current I A within 
a time interval At indicates with high probability that 
system B is occupied, and measuring a non-vanishing 
detector current indicates an empty channel B. 

To investigate within which limits on T, 77, and At 
this simplistic view is valid, we introduce the coarse- 
grained superoperators J LL = J (aQ \At), J LH = 
E„ B >i^ ( °" s) (At), J HL = E nA >iJ {nA0) m, and 
jHH - ^ nB>i l /(™>i»B)(At), corresponding to a van- 
ishing (L) or non-vanishing (H) current through channels 
A and B. That is, the L/iJ-discrimination threshold is 
set here independent of At. From Eq. (J2J), we may an- 
alytically obtain their Laplace transforms by performing 
the summation in frequency space (not shown) . We may 
then create coarse-grained current trajectories with both 
detectors A and B (leading to four different measurement 
outcomes 9Jt £ {LL, LH, HL, HH}), see inset of Fig. H 
Obviously, the correlation between blips of current I A 



and spikes of current I B is not perfect: The detector re- 
sult may fail, as e. g. tunneling charges through channel 
B may be missed. 

Without the FCS of channel B, the error probability 
of the detector can be calculated from the averaged joint 
probabilities of measuring a large current in channel A 
and a particle in channel B and measuring a vanishing 
current through A in combination with no particle in B: 
P e i(Af) = Tr{d B d B j£p} + Tr{(l-d B d B )J%p}, where 
p is defined via £{0,0)p = 0, Jj* = J HL + J HH , and 

J a = J r ' 



LL 



J . Its Laplace transform equates to 
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with v num (z) = 277r 2 (77 + r) + 2r(?7 + r)(4?7 + r)* 

+ (277 + r)(277 + 5r)z 2 +4(77 + r)z 3 + z 4 . From the Laplace 
transform we can directly conclude that for too short or 
too long measurement times the detector result will be 
useless: limAt-^o P<£r(At) = \ and liniAt->.oo Perr(At) = 
i . An intuitive explanation This can be understood intu- 
itively: For too short measurement times, the two peaks 
of the corresponding bimodal distribution will not have 
well separated, and measuring a vanishing particle num- 
ber (current) n A (At) is still possible even when channel 
B is empty. In contrast, for too large measurement times, 
we may obtain an average current over several cycles of 
loaded and empty channel B. For intermediate measure- 
ment times At however, we observe a pronounced mini- 
mum of P£ T {At), see Fig. [3] 

Dynamical Channel Blockade (Fig. Wf ). Here, 
both probed system (B) and detector (A) may be seen 
as bistable systems^, depending on whether 77 or T 
dominate, respectively. When occupied, both trans- 
port channels block the current through the other one 
completely^. With the replacements ft {(-a) = 1, 
ff(e A ) = 0, ff(e B ) = 1, and /£ (e B ) = in Eq. ©, it 
becomes obvious that the doubly charged state |11) will 
not be occupied. When we choose the junctions to A 
and B and the corresponding tunneling rates identically 
(77 —> r and £ — > %), we recover the formal Liouvillian 
structure in Ref.— for a two-terminal two-level system. 
The Laplace transform of the MGF M(x, £, z) may also 
be obtained analytically, and we find that the interac- 
tions and back-actions between system and detector are 
not negligible, exemplified by M(0, £, z) 7^ M B (£,z) and 
M (x, 0, z) ^ M A (x, z). From the MGF we obtain for the 
first cumulants {{n A )) = ^, {{n B )) = ^, and the second 
cumulants approach the asymptotic limits 
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such that the cross-correlations do now also appear in the 
currents. We may also here introduce the superoperators 
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FIG. 3: (Color Online) Detector error probability in the quan- 
tum detector limit [Eq. ©, solid lines] and in the dynamical 
channel blockade limit [Eq. JS), dashed lines] for varying ra- 
tios of average tunneled detector and system charges n/T. For 
a too small or too large number of tunneled detector charges, 
the detector result is practically useless, whereas for an inter- 
mediate number, a pronounced minimum (note the logarith- 
mic plot) of the error probability is found. Thin curves in the 
background mark the trajectory of the minimum. Asymptotic 
formula refer to solid lines (c) - dashed curve asymptotics 
would have slightly different prefactors. Filled symbols yield 
rough error estimates for experiments. 

corresponding to either zero or more than zero particles 
leaving either transport channel, and the current trajec- 
tories (not shown) are similar to the inset of Fig. [2j The 
Laplace transform P^ r (Ai) is obtained similarly 

pa ( v z 2 (3r + 2z) + 7 ? (r 2 + 6rz + 4z 2 ) 

eTl{ } ' 3z(T + z)[z(T + z)+r 1 (T + 2z)} 1 [ ! 

and we can directly conclude that for too short or too 
long measurement times the detector result will be useless 
since lim At ^ P c f r (At) = § and lim A t-!.oo P£ T (At) = §. 
For finite measurement times the error probability will 



show a pronounced minimum, sec Fig.[3j Therefore, mea- 
surements on A may give reliable information on the state 
of B (which differs strongly from its uncoupled dynamics) 

Detector Errors and Experiments. Eqns. ((SJ) 
and © generate a map of detector fidelity (Fig. [3J . The 
detection error only depends on the two dimensionless 
variables TAt (number of detector charges during the 
measurement time) and rj/T (ratio of system to detec- 
tor currents) and can be made extremely small. The 
striking similarity of the detector error probability for 
ideal and non-ideal detectors (also the similarity of the 
current trajectory in Fig. [2] with experiments for QPCs 8 ) 
indicates that the simple dependence on detector and sys- 
tem charge throughput is generic for 7™" <C P^t*- An 
error estimate taking into account a realistic At-scaling of 
the iiow/Jhigh discrimination threshold, asymmetric and 
energy-dependent tunneling rates, multiple levels etc. 
would be much more accurate, but for a crude estimate 
we may link experimental parameters with those of our 
model: The width of the experimental current divided 
by its mean value will in the high-bias limit yield the 
average number of tunneled detector charges during the 
measurement time via A/^V^St* ~ 1/VrAi, and the 
ratio of mean system current (number of blips divided by 
time) and mean detector current yields Pys/Pj 1 ^* ~ 
and we obtain error probabilities of approximately 10 -5 
per measurement for Refs^ ' 15 i 16 . 

Conclusions. Even with strong physical back-action 
the detector may yield reliable results for the instanta- 
neous occupation of channel B (with an altered dynam- 
ics) . To decide whether the used detector strongly back- 
acts on the probed system, it is therefore more useful 
to observe the scaling of the channel cross correlations: 
The smaller the scaling coefficient of the cross-cumulant 
in time, the smaller is the detector back-action. In order 
to estimate the fidelity of a single current measurement 
from Fig. [31 we propose to use the characteristics of time- 
resolved experimental detector current trajectories. 
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